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Abstract 

Recent advancement on the knowledge of multifragmentation and 
phase transition for hot nuclei is reported. It concerns i) the influence 
of radial collective energy on fragment partitions and the derivation of 
general properties of partitions in presence of such a collective energy, 
ii) a better knowledge of freeze-out properties obtained by means of a 
simulation based on all the available experimental information and iii) 
the quantitative study of the bimodal behaviour of the heaviest frag- 
ment charge distribution for fragmenting hot heavy quasi-projectiles 
which allows, for the first time, to estimate the latent heat of the 
phase transition. 
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1 Introduction 



Nucleus-nucleus collisions at intermediate energies offer various possibilities 
to produce hot nuclei which undergo a break-up into smaller pieces, which is 
called multifragmentation. The measured fragment properties are expected 
to reveal the existence on a phase transition for hot nuclei which was ear- 
lier theoretically predicted for nuclear matter [H EJ [3]. By comparing in 
detail the properties of fragments (Z > 5) emitted by hot nuclei formed 
in central (quasi-fused systems, QF, from 129 Xe+ nat Sn, 25-50 AMeV) and 
semi-peripheral collisions (quasi-projectiles, QP, from 197 Au+ 197 Au, 80 and 
100 AMeV), i.e. with different dynamical conditions for their formation, the 
role of radial collective energy on partitions is emphasized [1] and general 
properties of partitions are deduced (section 2). Then, in section 3, freeze- 
out properties of multifragmentation events produced in central collisions 
( 129 Xe+ nat Sn) are estimated [5] and confirm the existence of a limiting ex- 
citation energy for fragments around 3.0-3.5 MeV per nucleon. The deduced 
freeze-out volumes are used as a calibration to calculate freeze-out volumes 
for QP sources; thus one can locate where the different sources break in the 
phase diagram. Finally, in section 4, the charge distribution of the heaviest 
fragment detected in the decay of QP sources is observed to be bimodal. 
This feature is expected as a generic signal of phase transition in nonexten- 
sive systems such as finite systems. For the first time an estimate of the 
latent heat of the transition is also extracted [6] . 



2 Radial collective energy and fragment partitions 

To make a meaningful comparison of fragment properties which can be re- 
lated to the phase diagram, hot nuclei showing, to a certain extent, statis- 
tical emission features must be selected. For central collisions (QF events) 
one selects complete and compact events in velocity space (constraint of 
flow angle > 60°). For peripheral collisions (QP subevents) the selection 
method applied to quasi-projectiles minimizes the contribution of dynami- 
cal emissions by imposing a compacity of fragments in velocity space. The 
excitation energies of the different hot nuclei produced are calculated using 
the calorimetry procedure (see [3] for details). By comparing the properties 
of selected sources on the same excitation energy domain, significant differ- 
ences are observed above 5 AMeV excitation energy. QF sources have larger 
mean fragment multiplicities, < Mf rag >, even normalized to the sizes of 
the sources (which differ by about 20% for QF and QP sources), and lower 
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Figure 1: Evolution of the av- 
erage fragment multiplicity normal- 
ized to the source charge/size < 

>= < Mf rag JZ s > as a func- 
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tion of the relative velocity of frag- 
ments, p^el' ( see text) for different 
total excitation energy per nucleon 
of the sources. Full squares, open and 
full circles stand respectively for QF 
sources and QP sources produced at 
80 and 100 McV/nuclcon incident en- 
ergies. 



values for generalized asymmetry: Az = a z / ({Z) \/ M f rag — 1). A possible 
explanation of those different fragment partitions is related to the different 
dynamical constraints applied to the hot nuclei produced: a compression- 
expansion cycle for central collisions and a more gentle friction-abrasion 
process for peripheral ones. 
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Figure 2: Left and right sides re- 
fer respectively to the mean charge 
of the heaviest fragment of parti- 
tions, < Z\ >, and to the general- 
ized asymmetry in charge of the frag- 
ment partitions without the heaviest 
one, Az\{Zi}, (see text) as a func- 
tion of the reduced fragment muti- 
plicity, AIf rag /< Z s >, for different 
total excitation energy per nucleon 
of the sources. Full squares, open and 
full circles stand respectively for QF 
sources and QP sources produced at 
80 and 100 MeV/nucleon incident en- 
ergies. 
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The occurence of radial collective energy following a compression phase 
is predicted in semi-classical simulations of central collisions in the Fermi 
energy domain |8]. In experiments it was recognized, in most of the 



cases, from comparisons of kinetic properties of fragments with models. The 
mean relative velocity between fragments, (3 re i, independent of the reference 
frame, allows to compare radial collective energy for both types of sources 
(QF or QP). The effect of the source size (Coulomb contribution on frag- 
ment velocities) can be removed by using a simple normalization which takes 
into account, event by event, the Coulomb influence, in velocity space, of 
the mean fragment charge, {Z), on the complement of the source charge 

(Z s — (Z)); = Prei/ \/ {Z)(Z S — (Z)). At an excitation energy of about 

5 AMeV, the f}^2 values corresponding to QF and QP sources are similar. 
Above that excitation energy, values for QF sources exhibit a strong lin- 
ear increase, whreras for QP sources p^) only slightly increases up to 9-10 
AMeV. That fast divergence between the values of p£2 for the two types of 
sources signals the well known onset of radial collective expansion for central 
collisions. In [9j, estimates of radial collective energy (from 0.5 to 2.2 AMeV) 
for QF sources produced by Xe+Sn collisions are reported for four incident 
energies: 32, 39, 45 and 50 AMeV. Those estimates which were extracted 
from comparisons with the statistical model SMM assuming a self similar 
expansion energy have been used to calibrate the /3~f[ observable (see [I] for 
details). Then, one can show that radial collective energy is essentially due 
to thermal pressure for QP sources in semi-peripheral heavy-ion collisions as 
it is in hadron-induced reactions [10J. For QF sources produced in central 
heavy-ion collisions the contribution from the compression-expansion cycle 
becomes more and more important as the incident energy increases. Fig- 
ure [D shows, for different total excitation energy per nucleon corresponding 
to values defined by QF sources, the evolution of the average fragment multi- 
plicity normalized to the source charge/size < >= < Mf rag /Z s > as a 
function of f3^ el . We observe a well defined correlation which fully confirms 
the role of collective energy in producing more fragments. Depending on 
the source type the relative contributions to radial collective energy of ther- 
mal pressure and cycle compression-expansion strongly differ [3] but anyhow 
f5^ el , representative of the total collective energy, fixes the average degree of 
fragmentation (normalized mean fragment multiplicities). 

Does the intensity of the radial collective energy also govern the de- 
tails of fragment partitions, namely the relative charge/size of fragments 
in partitions. One can first consider the evolution of the size of the heavi- 
est fragment, for given total excitation energies, with the reduced fragment 
multiplicities Mf rag /< Z s >. On the left panel of fig. [2] average values of 
the heaviest fragment charge for QP and QF sources are reported: they 
follow exactly the same evolution. Finally the division of the charge among 



other fragments is investigated using the generalized asymmetry in charge 
of the fragment partitions. One can re-calculate the generalized asymmetry 
by removing Z\ from partitions, noted Az\{Z{\. The results, displayed in 
fig. [2] (right panel) do not depend on the source type at a given total ex- 
citation energy and a given reduced fragment multiplicity. Note that the 
general asymmetry follows a linear trend except the lower reduced fragment 
multiplicity which corresponds to Mf rag =2; indeed in that case, after re- 
moving Z\, only one fragment is available for the asymmetry calculation 
and in each event the heaviest Z of the partition below Z equal 5 was taken. 
Such a result shows the subtle role played by the radial collective energy. 
It influences the overall degree of fragmentation but it does not affect the 
relative size of fragments in partitions for fixed reduced fragment multiplic- 
ities. Those results represent a benchmark against which models describing 
fragmentation of finite systems should be tested. 

3 Freeze-out properties 

Starting from all the available experimental information of selected QF 
sources produced in central 129 Xe+ nai Sn collisions which undergo multifrag- 
mentation, a simulation was performed to reconstruct freeze-out properties 
event by event [5]. The method requires data with a very high degree 
of completeness, which is crucial for a good estimate of Coulomb energy. 
The parameters of the simulation were fixed in a consistent way includ- 
ing experimental partitions, kinetic properties and the related calorimetry. 
The necessity of introducing a limiting temperature related to the vanishing 
of level density for fragments |12j in the simulation was confirmed for all 
incident energies. This naturally leads to a limitation of their excitation 
energy around 3.0-3.5 AMeV as observed in [13 1 . The agreement between 
experimental and simulated velocity spectra for fragments of given charges 
(Z=6, 11, 18 and 27), for the different beam energies, is quite remarkable 
(see [5]). Finally relative velocities between fragment pairs were also com- 
pared through reduced relative velocity correlation functions [HI [151 EH E! 
(see fig. [3]). In the simulation the fragment emission time is by definition 
equal to zero and correlation functions are consequently only sensitive to the 
spatial arrangement of fragments at break-up and the radial collective energy 
involved (hole at low reduced relative velocity), to source sizes/charges and 
to excitation energy of the sources (more or less pronounced bump at v re d= 
0.02-0. 03c). Again a reasonable agreement is obtained between experimen- 
tal data and simulations, especially at 39 and 45 AMeV incident energies, 
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Figure 3: Comparison between the experimental (full points) and simulated (his- 
tograms) reduced relative velocity correlation functions for all the fragments. The 
reduced relative velocity between two fragments with charges Zi and Zj (Z^j >4) 
is defined as v re d—v re i/(Zi + Zj) 1 / 2 . Each panel refers to a different beam energy: 
32 AMeV (top left), 39 AMeV (top right), 45 AMeV (bottom left) and 50 AMeV 
(bottom right). From [5]. 



which indicates that the retained method and parameters are sufficiently 
relevant to correctly describe freeze-out topologies and properties. 

The major properties of the freeze-out configurations thus derived are 
the following: an important increase, from ~20% to ~60%, of the percentage 
of particles present at freeze-out between 32 and 45-50 AMeV incident en- 
ergies accompanied by a weak increase of the freeze-out volume which tends 
to saturate at high excitation energy. Finally, to check the overall physical 
coherence of the developed approach, a detailed comparison with a micro- 
canonical statistical model (MMM) was done. The degree of agreement, 
which was found acceptable, confirms the main results and gives confidence 
in using those reconstructed freeze-out events for further studies as it is done 
in [3]. Estimates of freeze-out volumes for QF sources produced in Xe+Sn 
collisions for incident energies between 32 and 50 AMeV evolve from 3.9 to 
5.7 V/Vq, where Vo would correspond to the volume of the source at normal 
density [5]. 



To calibrate the freeze-out volumes for other sources, we use the charge 
of the heaviest fragment < Z^ N ' > or the fragment multiplicity < >, 
normalized to the size of the source, as representative of the volume or 
density at break-up. From the four points for QF sources and the additional 
constraint that z[ N ^ = Mf rag =l at V/Vq=1, we obtain two relations V/Vq = 

fi(z[ N ^) and V/Vq = f2{Mj^J ), from which we calculate the volumes for 
QF sources at 25 AMeV and for QP sources. The results are plotted in fig-IU 
with error bars coming from the difference between the two estimates using 
/i and f2\ note that error bars for the QP volumes are small up to 7 AMeV, 
and can not be estimated above, due to the fall of < > at high energy 

(see fig. 5 of |4J). So only < (z[ >) can be used over the whole excitation 
energy range considered and the derived function is the following: 
V/V = exp(2A7 - 4.47 < {z[ N) >) + 0.86. 

The volumes of QP sources are smaller than those of QF sources (by 
about 20% on the E* range 5-10 AMeV). 
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Figure 4: Fragmentation position in 
the excitation energy-freeze-out vol- 
ume plane. The four full squares (QF 
sources) are taken from [5] . The open 
rectangle gives the estimated posi- 
tion (with error bar) for QF source at 
25 AMeV, and the open circles those 
for QP sources. From [4]. 



Z\ also presents some specific dynamical properties. As shown in [181117] 
for QF sources, its average kinetic energy is smaller than that of other 
fragments with the same charge. The effect was observed whatever the 
fragment multiplicity for Xe+Sn between 32 and 50 AMeV and for Gd+U 
at 36 AMeV. The fragment-fragment correlation functions are also different 
when one of the two fragments is Z\. This observation was connected to 
the event topology at freeze-out, the heavier fragments being systematically 
closer to the centre of mass than the others. 



4 Bimodality of the heaviest fragment and latent 
heat of the transition 



At a first-order phase transition, the distribution of the order parameter in 
a finite system presents a characteristic bimodal behaviour in the canonical 
or grandcanonical ensemble [19j . The bimodality comes from an anomalous 
convexity of the underlying microcanonical entropy [20] . It physically corre- 
sponds to the simultaneous presence of two different classes of physical states 
for the same value of the control parameter, and can survive at the ther- 
modynamic limit in a large class of physical systems subject to long-range 
interactions [21]. In the case of hot nuclei which undergo multifragmenta- 



Figure 5: Size of the heaviest frag- 
ment versus total excitation energy 
in AMeV. That picture is constructed 
using the fit parameters extracted 
from the equivalent-canonical distri- 
bution. The distance between the two 
maxima, liquid and gas peaks, pro- 
jected on the excitation energy axis 
corresponds to the latent heat of the 
transition. 

tion, the size/charge of the heaviest fragment was early recognized as an 
order parameter [22\ [23] using the universal fluctuation theory. In a quan- 
titative analysis for QP sources, the robustness of the signal of bimodality 
is tested against two different QP selection methods [6] . A weighting proce- 
dure p4j is used to test the independence of the decay from the dynamics of 
the entrance channel and to allow a comparison with canonical expectations. 
Finally, a double saddle-point approximation is applied to extract from the 
measured data an equivalent-canonical distribution. To take into account 
the small variations of the source size, the charge of the heaviest fragment 
Z\ has been normalized to the source size. After the weighting procedure, a 
bimodal behaviour of the largest fragment charge distribution is observed for 
both selection methods. Those weighted experimental distributions can be 
fitted with an analytic function (see [6] for more details) . From the obtained 
parameter values one can estimate the latent heat of the transition of the hot 
heavy nuclei studied (Z~70) as AE = 8.1(±0.4) atat (+1.2 - 0.9) syst AMeV. 
Statistical error was derived from experimental statistics and systematic er- 
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rors from the comparison between the different QP selections. The results 
(for one QP source selection) are illustrated in fig. [5j A detailed presentation 
and discussion of those results is also found in [25] . 

5 Conclusion 

Today a rather coherent and complete picture has been reached for a few 
exhaustive studies concerning multifragmentation and the related liquid- 
gas type phase transition which occurs at excitation energies between 2-3 
and 9-10 AMeV. Only the mechanism of fragment formation is still an open 
question: spinodal fluctuations in stochastic mean field approaches or many- 
body correlations early built in molecular dynamics [3]. At present, with 
the introduction of the N/Z degree of freedom in such studies, new signals 
predicted by theory must be investigated to precise and strengthen our ac- 
tual knowledge. It concerns fractionation (liquid more symmetric due to 
minimization of symmetry energy in dense phase), increased fractionation 
if spinodal instabilities are responsible for fragment formation, and the re- 
duction of spinodal zone for large N/Z values involved. For the future the 
introduction of the N/Z degree of freedom will permit to improve funda- 
mental information on the phase diagram of baryonic matter. Moreover hot 
exotic nuclei appear as a unique laboratory to serve as test-bench for theory 
of phase transitions of finite quantum systems with two components. 
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